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I.  SUMMARY 

A  three  year  study  of  the  effect  of  initial  conditions  on  the  topology  and  fine  scale 
structure  of  temporally  evolving  wakes  has  been  completed.  Early  in  the  project  new 
topological  analysis  tools  were  developed  and  used  to  analyze  available  simulations  of  the 
plane  mixing  layer.  During  the  course  of  these  preliminary  studies  it  was  found  that 
certain  geometrical  features  of  the  fine  scale  motions  in  the  mixing  layer  were  very  similar 
to  those  seen  in  simulations  of  forced  isotropic  turbulence.  As  a  consequence,  we 
adopted  a  two-track  approach  to  the  work;  one  track  was  to  gather  together  simulation 
data  on  as  wide  a  variety  of  shear  flows  as  possible  with  the  objective  of  identifying 
universal  feature  of  the  fine  scale  structure.  The  second  track  was  to  carry  out  a  new 
direct  simulation  of  the  plane  wake  in  which  systematic  vai'iations  of  the  initial  conditions 
were  used  to  identify  structural  features  of  the  developed  flow  which  depend  on  initial 
conditions. 

The  fundamental  questions  we  tried  to  address  were:  What  features  of  the  fine  scale 
motions  in  turbulent  flow  are  universal,  what  features  are  not,  and,  in  the  case  of  the 
plane  wake,  what  information  about  the  near  wake  survives  far  downstream? 
Inteipretation  of  the  computed  flow  fields  was  carried  out  using  topological  methods.  In 
this  approach  the  time  evolution  of  the  flow  is  followed  in  terms  of  joint  pdf’s  of  the 
invariants  of  the  velocity  gradient  ten,sor.  The  invariants  can  be  directly  related  to  the 
topology  of  the  local  velocity  field  enabling  useful  comparisons  to  be  made  between 
significant  features  of  the  pdfs  and  corresponding  features  in  physical  space.  Of 
particular  interest  is  the  second  invariant  of  the  rate-of-strain  tensor  which  is  directly 
proportional  to  the  dissipation  of  kinetic  energy.  An  intemediate  asymptotic  theory  for 
the  evolution  of  the  velocity  gradient  tensor  was  developed  under  this  project  and  u.sed  to 
identify,  for  the  first  time,  an  attractor  in  the  space  of  tensor  invariants.  Several  of  the 
important  universal  as  well  as  initial-condition-dependent  features  of  the  invariant  joint 
pdfs  were  explained  by  this  theory. 

II.  LONG  TERM  GOALS: 

1 .  To  gain  improved  understanding  of  the  fundamental  mechanisms  of  kinetic  energy 
dissipation  in  turbulent  flow  and  to  identify  universal  features  of  turbulent  fine  scale 
structure. 

2.  To  develop  dynamical  models  for  the  time  evolution  of  the  velocity  gradient  tensor. 

?>.  To  study  the  linkage  between  the  late  stages  of  flow  development  and  initial 

conditions. 

4.  To  develop  topological  methods  for  revealing  significant  flow  field  features 
contained  in  large  data  .sets 
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III.  APPROACH,  TASKS  COMPLETED 

I .  Method  of  approach,  sui-vev  of  simulation  data 

Studies  of  the  topology  of  fine  scale  motions  in  a  variety  of  Hows  have  been  can  ied 
out  during  this  project.  The  ca.ses  of  interest  include  the  following  time  developing 
simulations 


(i)  Weakly  compressible  mixing  layer 

(ii)  Incompressible  mixing  layer 

(iii)  Weakly  compressible  wake 

(iv)  Incompressible  wake 

(v)  Weakly  compressible  homogeneous  shear  fiow 

(vi)  Incompressible  homogeneous  shear  fiow 

(vii)  Incompres.sible  channel  fiow 

Discussion  of  these  cases  may  be  found  in  the  list  of  publications  attached  to  this  report 
(References  4,  5,  7.  8,  11,  and  12).  In  the  approach  the  partial  derivatives  of  the  velocity 
field  are  determined  at  every  point  in  the  flow.  These  are  u.sed  to  construct  the  invariants 
of  the  velocity  gradient  tensor  (P.  Q,  R),  the  rate-of-strain  tensor  (Ps,  Qs.  Rs)  and  the 
rate  of  rotation  tensor  (0,  Qw,  0).  For  incompressible  flow  the  first  invariant  is  zero. 
For  the  conditions  of  the  compressible  simulations,  the  first  invariant  is  found  to  be 
evei^where  small  relative  to  the  second  and  third  invariants  and  .so  in  both  cases  the  local 
topology  at  a  point  is  mainly  deteirnined  by  the  second  and  third  invariants.  The  resulting 
en.semble  of  (Q,R),  (Qs,  Rs)  and  (Qs,  Qw)  pairs  is  sorted  into  bins  and  presented  as 
number  density  contours.  See  for  example  Figures  7  to  18  in  Reference  4.  In  this  forni 
the  invariant  plots  are  equivalent  to  unnormalized  joint  probability  density  disuibutions  of 
the  two  invariants  being  crossplotted.  The  invariants  can  be  directly  related  to  the 
topology  of  the  local  velocity  field  enabling  useful  comparisons  to  be  made  between 
significant  features  of  the  joint  pdfs  and  corresponding  flow  structures  in  physical  space. 
Working  in  the  space  of  invariants  in  this  fashion  is  an  efficient  way  of  revealing 
significant  structural  features  of  the  flow.  Of  particular  interest  is  the  second  invariant  of 
the  rate-of-strain  tensor,  Qs  ,  which  is  directly  proportional  to  the  dissipation  of  kinetic 
energy. 

2.  Modeling  the  evolution  of  the  veKx:itv  gradient  ten.sor 

The  velocity  gradient  tensor  satisfies  a  nonlinear  evolution  equation  of  the  forni 
dAy/dt  -I-  Aj^Ajjj  -  (1/3)  where  A.j=3uj/3xj  and  the  ten.sor  Hy  contains 

terms  involving  the  action  of  cross  derivatives  of  the  pressure  field  and  viscous  diffusion 
of  the  velocity  gradient.  The  homogeneous  case  (Hy  =  0)  forms  a  ninth  order  system  of 
ODE'S  for  the  components  of  the  velocity  gradient  ten.sor  and  has  been  studied  by  a 
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number  of  investigators  using  asymptotic  analysis  and  numerical  integration  beginning  in 
the  early  198()'s.  In  the  course  of  our  studies  we  discovered  an  exact  integral  of  the 
equations  (Reference  1).  The  asymptotic  behavior  of  the  exact  solution  has  a  number  of 
properties  in  common  with  observations  of  direct  simulation.  In  addition  the  exact 
solution  has  enabled  us  to  precisely  determine  the  connection  between  initial  conditions 
and  the  asymptotic  state  of  the  restricted  Euler  system. 

No  one  previously  has  been  able  to  make  any  general  statements  about  the  behavior 
of  the  velocity  gradient  ten.sor  when  the  as.sumption  H.j  =  0  is  removed.  However  using 
re.sulis  from  the  restricted  Euler  solution  a  new  intermediate  asymptotic  theory  for  the 
evolution  of  Ajj  has  been  developed  in  which  no  restriction  is  placed  on  but  the 
velocity  gradient  tensor  is  assumed  to  evolve  toward  the  form  dAjj/dt=f(t)Ajj  where  f(t)  is 
a  scalar  function  of  time.  The  restricted  Euler  solution  predicts  a  specific  f(t)  however  in 
the  context  of  the  intennediate  a.symptotic  model  f(t)  is  unspecified  and  can  be  used  to  fit 
experimental  (numerical)  data.  Under  this  assumption  a  general  relationship  can  be 
worked  out  between  the  invariants  Q  and  R  and  the  discriminant  of  H-j .  This  leads  to  the 
notion  that  Q  and  R  must  lie  in  a  region  of  attraction  in  the  space  of  tensor  invariants  (See 
Figures  ?>  and  4  of  Reference  3).  Various  regions  of  this  attractor  represent  local  fiow 
fields  which  can  be  characteri/ed  as  either  sheet-iike  or  tube  like  in  nature. 

3.  Direct  simulation  of  temp(>rallv  evolving  planar  incompressible  wakes 

Previous  experimental  studies  of  plane  wakes  produced  by  an  assortment  of 
different  wake  generators,  all  arranged  to  have  the  same  drag,  demonstrated  strong 
.sensitivity  of  the  turbulent  .structure  of  the  wake  to  initial  conditions.  This  sensitivity  is 
observed,  to  varying  degrees,  in  a  wide  variety  of  free  shear  Hows  and  is  one  of  the  most 
difficult  impedimenLs  to  the  development  of  accurate  models.  In  the  pre.sent  study  direct 
numerical  simulations  of  an  incomprc.s.sible  plane  wake  are  used  to  investigate  the  linkage 
between  the  late  stages  of  wake  development  and  properties  of  the  initial  fiow  field.  A 
.series  of  simulations  were  carried  out  on  the  NAS  Cray  YMP  to  explore  the  effects  of 
various  combinations  and  phasings  of  2-D  and  3-D  laminar  disturbances  on  the  Uansition 
to  turbulence  in  wakes  with  Reynolds  numbers  based  on  initial  halfwidth  and  deficit 
velocity  of  between  5(K)  and  4(KK).  The  simulations  used  a  pseudospectral  Fourier/Jacobi 
code  developed  at  NASA  Ames.  All  the  simulations  were  initiated  with  a  Gau.ssian  mean 
velocity  field  perturbed  by  small  disturbances  periodic  in  the  streamwise  and  spanwi.se 
directions.  Some  of  the  preliminary  results  of  these  studies  are  described  in  Reference 
10.  Scaling  laws  and  the  effects  of  initial  conditions  on  the  structure  of  invariant  joint 
pdfs  are  pre.sently  under  study. 
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4.  Development  of  a  direct  Fourier  method  for  simulating  incompressihle  Hows  with 

one  infinite  and  two  periodic  directions 

Initially  the  simulation  code  for  the  wake  was  a  modification  of  an  existing  code 
developed  for  the  Cray-YMP.  However  because  of  the  high  level  of  usage  of  the  YMP 
the  turn  around  time  for  the  computation  was  exuemely  slow  and  it  was  decided  to  move 
the  simulation  to  the  iPSC/86()  hypercube.  An  entirely  new  spectral  code  has  been 
written  using  a  novel  approach  in  which  a  Fourier  approximation  to  the  How  in  the 
vortical  region  is  matched  to  an  analytic  asymptotic  solution  for  the  far  field.  The 
algorithm  makes  use  of  the  fact  that  the  vorticity  field  in  planar  free  shear  Hows  is 
compact  in  the  cross-stream  direction.  This  allows  the  vorticity  equations  governing  the 
How  to  be  solved  in  a  fully  periodic  box,  with  the  non-periodic  velocity  being  calculated 
by  matching  to  tne  known  a.symptotic  solution.  Because  this  avoids  nonlinear  mapping 
of  the  How  domain,  direct  use  of  Fast  Fourier  Transforms  in  all  three  spatial  directions  is 
possible,  resulting  in  a  code  which  is  small,  fast,  and  natural  to  implement  on  parallel 
architectures.  Further  studies  of  the  plane  wake  using  this  code  along  with  comparisons 
with  the  previous  Cray  YMP  simulations  are  presently  in  progress. 

IV.  SCIENTIFIC  /  ENGINEERING  RESULTS 

The  work  carried  out  under  this  grant  has  revealed  a  number  of  important 
geometrical  features  of  the  fine  scale  structure  of  incompressible  and  weakly 
compressible  inhomogeneous  shear  Hows;  The  second  and  third  invariants  of  the 
velocity  gradient  and  rate-of-strain  tensors  are  strongly  correlated;  Motions 
corresponding  to  high  rates  of  dissipation  are  found  to  be  characterized  by  a  ."^-D  rate-of- 
strain  topology  which  is  of  the  type  saddle-saddle-unsiable-node;  Regions  conesponding 
to  high  rates  of  di.ssipation  are  found  to  be  characterized  by  comparable  magnitudes  of 
enstrophy  and  strain  with  the  vorticity  vector  closely  aligned  with  the  intermediate 
positive  strain  direction.  In  the  incompressible  cases,  motions  characterized  by  the 
highest  enstrophy  were  found  to  be  nearly  strain-free  (ie,  in  solid  body  rotation). 
Whereas  in  the  compressible  cases  motions  characterized  by  the  highest  enstrophy  were 
always  imbedded  in  regions  of  high  strain  rate.  By  studying  a  variety  of  Hows  it  has 
been  po.ssible  to  accumulate  a  large  body  of  information  which  supports  the  hypothesis 
that  in  free  shear  layers  the  features  itemized  above  are  universal  (References  4,  5,  6,  and 
7).  Recently  in  Reference  4  it  was  discovered  that  the  rate-of-strain  topology  just 
de.scribed  persists  to  substantially  larger  .scales  than  was  previously  realized. 

The  asymptotic  behavior  of  the  ivstricted  Euler  .solution  discu.s.sed  above  reproduces 
many  of  the  topological  features  seen  in  the  simulations.  Examination  of  the  exact 
.solution  of  this  system  provides  new  insight  into  the  manner  in  which  the  fine  scale 
motions  of  turbulence  can  evolve  to  a  universal  state  and  the  degree  to  which  this  state 
can  depend  on  initial  conditions.  A  new  inteiTnediale  asymptotic  model  of  the  evolution 
of  the  velocity  gradient  tensor  shows  improved  agreement  with  simulation  data  and 
identifies,  for  the  first  time,  a  region  of  attraction  in  the  space  of  velocity  gradient  tensor 
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invariants.  Very  recently  in  a  study  of  turbulent  channel  How  (Reference  5)  it  has  been 
found  that  in  regions  of  high  dissipation  the  voiticity  vector  is  aligned  approximately  10 
degrees  away  from  the  intermediate  strain  direction. 

Simulations  of  plane  incompressible  wakes  demonstrate  the  important  role  of 
subharmonic  components  in  the  initial  field  for  the  subsequent  generation  of  three- 
dimensionality. 

V.  ACCOMPLISHMENTS 

1 .  A  three-dimensional  topological  analysis  method  based  on  the  joint  pdf  of  local 
tensor  invar  iants  has  been  developed  and  found  to  be  a  useful  way  to  elucv'aie  significant 
features  of  large  fiow  field  data  sets. 


2.  The  topological  analysis  method  has  been  used  to  identify  important  geometrical 
features  of  the  fine  .scale  structure  of  turbulence  in  a  variety  of  incompressible  and  weakly 
compre.ssible  inhomogeneous  shear  Hows.  Evidence  to  date  supports  die  hypothesis  that 
the.se  features  are  universal. 


?>.  An  exact  clo.sed  form  solution  has  been  di.scovered  for  a  restricted  Euler  equation 
for  the  velocity  gradient  tensor  giving  new  theoretical  insight  into  the  evolution  of 
turbulent  fine  .scale  structure  to  a  univer.sal  state.  An  intermediate  asymptotic  theory 
suggests  for  the  first  time  that  the  velocity  gradient  tensor  evolves  to  an  attractor  in  the 
.space  of  tensor  invariants. 

4.  Simulations  of  plane  incompr'essible  wakes  demonsti'ate  the  importance  of 
subharmonic  components  in  the  initial  field  for  the  subsequent  generation  of  three- 
dimensionality.  A  new  spectral  code  has  been  developed  in  which  a  Fourier 
approximation  to  the  flow  in  the  vortical  region  is  matched  to  an  analytic  asymptotic 
solution  for  the  far  field.  This  allows  direct  use  of  FFT’s  in  all  three  spatial  directions, 
resulting  in  a  code  which  is  especially  well  .suited  to  parallel  architectures. 
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